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We report detailed inelastic neutron scattering measurements on single crystals of the frustrated 
two-leg ladder BiCu2P06, whose ground state is described as a spin liquid phase with no long-range 
order down to 6 K. Two branches of steeply dispersing long-lived spin excitations are observed 
with excitation gaps of Ai = 1.90(9) meV and A2 = 3.95(8) meV. Significant frustrating next- 
nearest neighbor interactions along the ladder leg drive the minimum of each excitation branch 
to incommensurate wavevectors £1 — 0.5747T and £2 = 0.5537T for the lower and upper energy 
branches respectively. The temperature dependence of the excitation spectrum near the gap energy 
is consistent with thermal activation into singly and doubly degenerate excited states. The observed 
magnetic excitation spectrum as well as earlier thermodynamic data could be consistently explained 
by the presence of strong anisotropic interactions in the ground state Hamiltonian. 

PACS numbers: 75.10.Jm,75.10.Kt,75.40.Gb 



I. INTRODUCTION 

Low-dimensional quantum antiferromagnets can real- 
ize a rich and diverse array of physical phenomena from 
a seemingly simple set of interactions. The canonical 
low-dimensional quantum antiferromagnet is the spin- 
1/2 Heisenberg chain, in which neighboring spins are 
coupled with antiferromagnetic exchange J\. The sys- 
tem does not have any long range order and the spin-spin 
correlations decay with a power law. In spin-1/2 chains 
the elementary excitations are S = 1/2 quasiparticles 
termed spinons and the dynamic susceptibility is dom- 
inated by a gapless dispersive continuum!^ Frustration 
can be introduced to the spin-1/2 chain by competing 
antiferromagnetic next-nearest-neighbor (NNN) interac- 
tion J 2 . For a large NNN exchange of J2/J1 > 0.241 the 
ground state is composed of dimerized sing and the 

excitation spectrum is described with a coherent triplet 
band of gapped excitations. At the Majumdar-Ghosh 
(MG) point J2/J1 = 0.5 the singlet ground state is exact 
and the system forms a one-dimensional dimer crystal. 5 
As the NNN exchange is increased beyond the MG point 
frustration drives the spin correlations in the dimerized 
system to an incommensurate wave vectorial 

A gap may also be introduced into the excitation spec- 
trum by coupling two neighboring chains via a short 
range interaction J rung , to create an even leg ladder 8 . 
The rung coupling confines spinons into S=l magnons 
and the low energy dynamic susceptibility exhibits a well 
denned, triply degenerate, single particle peakP Theoret- 
ical understanding of the properties of even-leg ladders 
is highly developecPE^] and a number of experimental re- 
alizations have allowed for a detailed understanding of 
the very rich physical phenomena present in spin-ladders 
for both the strong rung coupling (rung-singlet ^M^S and 



strong- leg coupling (Haldane) 15 * 161 regimes. However, the 
physics of spin ladders with additional frustrating NNN 
interactions is not understood well. Field theoretical and 
numerical investigations have indicated that the addition 
of frustration can lead to many exotic quantum ground 
states not realized in the standard two- leg ladder EZHED 
yet comprehensive experimental investigations are still 
lacking. 

Such a frustrated two-leg ladder seems to be real- 
ized in a newly discovered quas i- two-dimensional spin- 
1/2 compound BiC^POe! 20 1 21 1 The crystal structure 
of BiCu2P06 is shown in Fig. [I] Zig-zag chains of 
Cu 2+ ions run parallel to the crystallographic b-axis, as 
shown clearly by a projection in the a — b plane in Fig. [I] 
(c). Along the chains Cu 2+ ions at crystallographically 
inequivalent sites, labeled Cua and Cub, interact via 
NN antiferromagnetic exchange J\ and NNN antiferro- 
magnetic exchange J2. The chains are coupled along 
the c-axis by both J3 and J4. Although the J3 inter- 
action distance is shorter than J 4 , band structure calcu- 
lations predict J4 to be the dominant coupling along the 
c-axis. 23 Structurally, this is a result of the particular su- 
perexchange pathways. J3 is a nearly 90° Cu-O-Cu bond, 
which is usually small and ferromagnetic, in contrast J4 is 
mediated by a nearly 180° Cu-O-O-Cu bond. Therefore, 
BiCu2P06 can be considered as a system o f J\- J 2- J a lad- 
ders with a weaker interladder coupling J3! 21 * 23 * 

Magnetic susceptibility, magnetization, and heat ca- 
pacity measurements on powder samples have consis- 
tently reported that BiCu2P06 has a singlet groun d state 
with a spin excitation gap of ~2.9 meVpSESHH] The 
thermodynamic measurements cannot be interpreted in 
terms of simple gapped one-dimensional (ID) mod- 
els, including the J1-J2 chain and the two- leg ladder, 
but require a combination of frustration and two-leg- 
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FIG. 1. (a) Schematic representation of the crystal structure 
of BiCu2P06- The unit cell is orthorhombic, space group 
Pnma with a = 11.755 A, b = 5.16 A, c = 7.79 A at 6 kP 

(b) Perspective view of the ladder unit in the b-c plane and 

(c) projection into the a-c plane showing the zigzag chains of 
Cu 2+ ions. 

ladder geometry! 21 * 24 * Inelastic neutron scattering (INS) 
measurements on powder samples have reported that 
BiCu2P06 has a dispersive spin excitation spectrum with 
a gap of A = 2.9 meV and that significant frustration 
is required to describe the powder averaged structure 
factor! 2 ^ More recent high field magnetization measure- 
ments on single crystals have revealed that the critical 
field for the closing of the spin gap is dependent on the di- 
rection of the applied field indicative of strong anisotropic 
interactions in the magnetic Hamiltonian! 2 ^ 

In real materials (usually) small anisotropic interac- 
tions, additional to the Heisenberg exchange couplings, 
invariably exist. Such interactions are often a negligi- 
ble perturbation to the isotropic Heisenberg interaction, 
as in cuprate super conduct or s^^Zl, but can sometimes 
alter the magnetic properties in a fundamental way. As 
an example, for compounds that lack local invers ion sy m- 
metries anisotropic Dzyaloshinksky-Moriya (DMp^ in- 
teractions are permitted in the magnetic Hamiltonian. 
Depending on the particular crystallographic symmetries 
the DM interaction can induce a staggered fieldP^ In ID 
chains this staggered field results in an effective confine- 
ment potential between spinons and the appearance of 
incommensurate gapped modes on application of a mag- 
netic field. ESI The staggered field also has a drastic effect 
for spin ladders. In contrast to a spin ladder in a uni- 
form magnetic field, which transitions to a gapless phase 
above a critical field, the presence of a staggered field is 
predicted to lead to a quantum phase transition between 
two gapped phases above and below the critical field.^ 

The numerous exchange pathways in 
BiCu2P06 complicate interpretation of thermody- 
namic data and important details including interladder 
coupling and anisotropic interactions are often neglected 
in the analysis. In order to determine the microscopic 
spin Hamiltonian of a complex magnetic system with 
competing interactions and anisotropic exchange cou- 



plings, INS measurements using a single crystal sample 
are essential. 

We have conducted an extensive neutron scattering 
study of the magnetic excitation spectrum in BiCu2P06- 
Our data confirms that BiCu2P06 is appropriately de- 
scribed by weakly interacting two-leg ladders with incom- 
mensurate dynamic correlations driven by frustration. In 
contrast to the single triply degenerate excitation branch 
expected for an isotropic ladder, we observe two branches 
of steeply dispersing, long-lived, excitations. The excita- 
tion gap of each mode, directly probed by INS, was mea- 
sured to be Ai = 1.90(9) meV and A 2 = 3.95(8) meV, 
this differs significantly from the ~ 2.9 meV gap ex- 
tracted from thermodynamic measurements^ 2 ^ 2 ^ indi- 
cating that current models are inadequate to describe 
the ground state of BiCu2P06- Furthermore, the tem- 
perature dependence of the incommensurate modes are 
consistent with a description in terms of thermal activa- 
tion into singly and doubly degenerate modes. We argue 
that, in addition to frustration, strong anisotropic inter- 
actions are important for understanding the physics of 
this material. 



II. EXPERIMENTAL DETAILS 

Experiments were carried out on a 4.5 g single crystal 
sample grown using the traveling floating zone method. 
The sample mosaic was measured by neutron scatter- 
ing to be 0.2° at T = 6 K. Measurements were per- 
formed on the C5 DUALSPEC triple axis spectrometer 
at the Canadian Neutron Beam Centre at Chalk River 
Laboratories and on the HB1 triple axis spectrometer 
at HFIR. Both instruments employed a vertically focus- 
ing pyrolytic graphite (PG) monochromator and C5 was 
equipped with a flat graphite analyzer while HB1 uti- 
lized a fixed vertically focusing analyzer. On C5 the 
sample was mounted in the (0, k, I) scattering plane and 
the spectrometer was operated at a fixed final energy of 
14.56 meV. Experiments on HB1 were performed with 
the sample mounted in the (h, k, 2k) horizontal scatter- 
ing plane at a fixed final energy of 14.7 meV. Tempera- 
ture control was provided by a closed cycle cryostat. All 
data was corrected for higher-order wavelength neutrons 
in the incident beam monitorP 2 ^ Intensities were placed 
on an absolute scale by normalization with the integrated 
intensity of a transverse acoustic phonon measured near 
the (004) Bragg peak on the respective instrument. 

The crystallographic unit cell of BiCu2P06 is or- 
thorhombic, space group Pnma with a = 11.755 A, b 
= 5.16 A, c = 7.79 A. Along the zig-zag chains NN 
Cu 2+ ions in BiCu2P06 are separated by b/2 so that the 
momentum of magnetic excitations in the b* direction is 
indexed using k = q • b/2 = 7rk. 
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III. EXPERIMENTAL RESULTS 

A. Spin excitation spectra 

The momentum and energy dependence of spin exci- 
tations in BiCu2PC>6 were surveyed through a series of 
constant-Q scans. No evidence for elastic magnetic scat- 
tering at T = 6 K was found indicating the absence of 
static magnetic order in BiCu2P06- Representative scans 
taken in proximity of the spin gap at T=6 K are shown 
in Fig. |2j Throughout most of the Brillouin zone the 
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FIG. 2. Representative constant-Q scans. Data in panels 
(a)-(d) were obtained on C5 with the sample aligned in the 
(0, k, I) scattering plane and a collimation of 33'-48'-51'-144\ 
Data in panels (e) and (f) was obtained on HB1 with the 
sample aligned in the (h,k, 2k) scattering plan and a collima- 
tion of 48'-40'-40'-120'. All data was collected at T = 6 K. 
The solid black lines are the results of a global fit to the sin- 
gle mode approximation, gold and blue filled areas show the 
contribution from each mode including resolution effects. 

scattering intensity is dominated by two well defined and 
highly dispersive modes. The dispersion of these modes 
is not commensurate with the structural unit cell and 
the scattering intensity vanishes as the temperature is 
increased above 60 K, where a broad maxima in the mag- 
netic susceptibility was observed^, confirming the mag- 
netic origin. Extensive sampling of Q-values throughout 
reciprocal space enabled the construction of a map of 
the dynamic structure factor 5(q, uj) shown in Fig. [3] (a) 
-(c). 

Magnetic excitations in BiCu2P06 are highly 
anisotropic, with a bandwidth of ~12 meV in the 
b* direction and ~2 meV in the c* direction. Dispersion 
along the c* direction confirms the presence of a sizable 
interladder interactions. Any dispersion along the a* 
direction could not be resolved by our thermal triple 
axis measurements, indicating that BiCu2P06 should 
be regarded as a quasi-two-dimensional system. 

Within the first Brillouin zone, each branch of the ex- 



citation spectrum has two minima occurring at wavevec- 
tors ( a = 0.5 + ^ and ( a — 1.5 — 5 a . The distinct 
double-well structure is reminiscent of the dispersion in 
J1-J2 chains and the incommensurate minima is indica - 
tive of a strong competition between J\ and ^ JHZHZHU 

The excitation spectrum shown in Fig. [3ja) cannot be 
described using strong-coupling expansions for a two-leg 
ladder with competing interactions ! 19 * 23 * These perturba- 
tive expansions require that the rung coupling (J4) is 
much larger than the coupling along the chains (Ji^)- 
The large bandwidth of each branch compared with the 
gap energies implies that BiCu2P06 likely lies in an in- 
termediate coupling regime where J\ ~ J4, this is in qual- 
itative agreement with band-structure calculations P2J 

Lacking a microscopically motivated theory for 
BiCu2P06, the essential features of the lowest energy ex- 
citations were determined by fitting data in the vicinity 
of the dispersion minimum with a single mode approxi- 
mation (SMA) cross-section and an empirical dispersion 
relation: 



S(q,w) = 



2 • 2 
v„ sin 



(1) 

(2) 



is a mode dependent intensity pre-factor, v a param- 
eterizes the spin- wave velocity of each mode, and we set 
h = ks = 1 throughout this paper. Equation ([2| was 
convolved with the instrumental resolution function 33 34 
and globally fit to scans in the range 0.57T < k < 0.7tt 
and 1.37T < k < 1.57T, only A^ was allowed to vary be- 
tween scans. Incoherent background was modeled as a 
constant plus a Gaussian function centered on w = 0. 
The broad peak widths and Q-dependence of lineshapes 
visible in Fig. [2] (a) -(d) and Fig. [S] (a)-(b) 

are accounted 

for entirely by resolution effects and two steeply dispers- 
ing modes indicating that the excitations are long-lived. 
Additional intensity appears in (/i, /c, 2k) plane which is 
best accounted for by a damped harmonic oscillator, blue 
shaded area in Fig. 21(e) and (f)P^ Dispersion parameters 
determined from the low-energy global fit were: Ai = 
1.90(9) meV, A 2 = 3.95(8) meV, v t = 16.4(1.4) meV, 
v 2 =19.6(1.7) meV, 5 t = 0.074(2), 5 2 = 0.053(4). Scans 
simulated using these parameters are shown as solid lines 
in Fig. [2] and the resulting dispersion is plotted as solid 
and broken grey lines in Fig. [3] (a). 

The incommensurate wavevector of the lower branch 
Ci = 0.5747T can be used to estimate the ratio J2/J1 
in the two limits of strong rung coupling J 4 > J\ and 
uncoupled weakly interacting J1-J2 chains. For strong 
rung coupling the wavevector minimizing the dispersion 
can be estimated exactly^ 



c 



-Ji 

4J 2 



(3) 



giving J2/J1 = 1.08. In the opposite limit of iso- 
lated J1-J2 chains, density matrix renormalization group 
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FIG. 3. (a) - (c) Inelastic neutron scattering intensity map for T = 6K constructed by linear interpolation of a series of constant 
scans. All data have been corrected for the isotropic Cu 2+ form factor. Solid and dashed gray lines show the dispersion relation, 
Eq. 0, using the parameters obtained from the single mode approximation fit. Solid white symbols in (c) show the fitted peak 
positions of resolution limited mode, (d) Structure factor data obtained at T = 6 K, for equal-time spin correlations along the 
chain b* direction, calculated by numerically integrating constant q scans for uj > 0.5 meV. (e) - (f) Measured first moment as 
a function of wavevector along the rung c* direction and perpendicular to the planes of the ladders, along a*. Solid lines are 
fits as described in the text. 



(DMRG) results 7 can be used to estimate J2/J1 from ( 
giving J2/J1 « 0.88. Our INS data situates BiCu 2 P0 6 in 
an intermediate coupling regime; however, the two lim- 
its of isolated chains and strong rung coupling tightly 
constrain J 2 « J\. 

Information about the length scales associated with 
spin correlations can be extracted from the equal-time 
structure factor, shown in Fig. [3] (d). It was calculated 
by numerically integrating constant q scans between 0.5 
and 15 meV, all scans were corrected for the isotropic 
Cu 2+ form-factor prior to integration.^ Static correla- 
tions in BiCu2P06 are well described by a double peaked 
square root Lorentzian, characteristic of spin chains with 
competing NN and NNN interactions and exponentially 
decaying spin-spin correlation^ 

S(k) oc 1 = + 1 = , (4) 

where C = 0.5 + 5, C" = 1.5 - J and f = 1/k is the 
correlation length. Fitting Eq. Q to data in Fig. [3] (d) 
yields a spin-spin correlation length of £/ (b/2) = 8.0(6) 
and 6 = 0.074(5). 

The ~ 2 meV bandwidth along the c* direction implies 
that there is a significant coupling between the ladder 
units. However, it is not possible to distinguish either J3 
or J 4 as the inter /intra ladder coupling from the shape 
of the dispersion alone. More information is contained 
in the modulation of INS intensity along the c* and a* 
directions. In the absence of a microscopic model for 



the spin dynamics in BiCu2P06 the first moment sum 
rule can be utilized to determine the relative contribution 
of spin-pair corre lations across each bond to the ground 
state energy!*™*" For a Heisenberg Hamiltonian the first 
moment sum rule is written 

/oc 
ujS aa (q,oj)doj 
-CO 

= - - J n' ( S i S i' ) I 1 - cos (<1 • r ii' )] , (5) 

where (Sfsf,) is the spin-spin correlation across bond 
rjjf. Within the single mode approximation the dimer- 
interference term [1 — cos (q • r^y)] accounts entirely for 
the q-dependent intensity modulation of the inelastic 
scattering intensity along the rung direction for a sys- 
tem of isolated ladder jJHESI However, the sum rule 
in Eq. ([5| is strictly only valid for a centro-symmetric 
lattice obeying inversion symmetry, hence must be ap- 
plied with caution when analyzing the spectrum from 
BiCu2P06 where a significant DM anisotropy on each 
rung is potentially present. The first-moment was mea- 
sured in BiCu2P06 by numerically integrating constant- 
q scans along the a* and c* directions and globally fit to 
equation All scans were corrected for the isotropic 
Cu 2+ form-factor prior to integration!^ The best fit was 
obtained with J 3 (S S 3 ) = -12(5) meV, J 4 (S S 4 ) = - 
58(2) meV. Intensity modulation along the c* and a* 
directions is well described including only bonds J3 and 
J4 the results of the fit are shown as solid lines in Fig. [3] 
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(e) - (f). The intensity modulation is thus consistent 
with a two-leg ladder formed by antiferromagnetic rung 
coupling J4 and weaker interladder coupling J3. 

Our neutron scattering results confirm that the J\- 
J2-J4 two-leg ladder with interladder exchange J3 is 
an appropriate description of BiCu2P06 . However, 
the gap energy differs significantly from the value of 
A ~ 2.9 me V repor ted by previous thermodynamic 
measurements.^ 21 * 23 * This discrepancy is potentially re- 
solved by the presence of strong anisotropic interactions 
neglected in the thermodynamic analysis. Indeed, DM 
anisotropics are permitted by the crystallographic sym- 
metries of BiCu2P06 . These anisotropics can act to 
split the degeneracy of the lowest lying excitations in 
zero field, so that the thermodynamic properties are con- 
trolled by two gaps. It is notable, if perhaps coinciden- 
tal, that the 2.9 meV thermodynamic gap is in agree- 
ment with the average of the two gaps measured by INS, 
l/2(Ai + A 2 ) = 2.93(6) meV. We have also checked that 
the low-temperature magnetic specific heat can be fit 
using an effective ID model including the sum of con- 
tributions from both the low and high energy modes 
C m = 1/2 (#1 exp(-Ai/T) +# 2 exp(-A 2 /T)) with each 
gap fixed at the value determined from INS. Further- 
more, a zero-field splitting of the lowest lying excitation 
resulting from a staggered DM interaction is consistent 
with the temperature dependence data presented next. 
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FIG. 4. Temperature dependence of spin excitations near the 
incommensurate wavevector at (0,1.425,1), data collected on 
C5. (a) A series of representative constant-q scans, data has 
been offset for clarity, solid lines are fit to resolution convolved 
cross-section, (b) - (d) Results of fit to resolution convolved 
Lorentzian scattering cross-section. Solid and dashed lines 
in (b) are fit to the population thermalization factor from 
RPA theory. The solid and dashed lines in (c) are a fit to 
exponential activated behaviour as described in the text. 



B. Temperature Dependence 

The temperature dependence of spin excitations at 
q = (0, 1.425, 1) near the incommensurate wavevector is 
shown in Fig. ma). As the temperature is increased the 
inelastic features broaden and shift to higher energies. 
Above 25 K the two excitations can no longer be dis- 
tinctly identified and above 55 K the inelastic signal can- 
not be distinguished from background. Each spectrum 
was fit to the resolution convolved single mode cross sec- 
tion [Eq. ([2])] with a normalized Lorentzian function in 
place of the delta function to account for intrinsic broad- 
ening of each mode. Lacking detailed measurements of 
the dispersion for T > 6 K the spin wave velocity was 
assumed to be independent of temperature for all fits. 
Constant-q scans were performed at a wavevector where 
the orientation of the resolution ellipsoid with respect 
to the slope of the dispersion minimizes any resolution 
effects. We have also checked that allowing for a tem- 
perature dependent spin wave velocity does not affect 
the results presented here. Figures [4] (b) - gc) show the 
temperature dependence of the gap energy, half-width 
half-maximum and integrated intensity of each mode. 

In accordance with the thermalization of the ground 
and excited states, the integrated intensity decreases 
with increasing temperature. The random phase ap- 
proximation (RPA) for weakly interacting dimers pre- 
dicts the scattering intensity to scale with the ground- 
state excited- state population difference An(A a /T) = 



(l-exp(-A a /T))/(l + 77exp(-A a /T)), where rj is 
the degeneracy of the level!™™! The RPA form is a good 
description of the data after fitting only an overall scale 
factor [solid lines in Fig. [4] (b)] . The best fit is obtained 
assuming the two modes are composed of a lower energy 
doubly degenerate excitation (Ai) and a higher energy 
singly degenerate excitation (A 2 ). Although we do not 
have sufficiently detailed data to completely rule out the 
possibility of two triply degenerate modes. The scaling 
of the intensities provides good qualitative evidence for 
the splitting of a single triply degenerate excitation band 
into a doublet and singly degenerate excitation. 

For systems with a gapped excitation spectrum in- 
creasing temperature results in an increased density 
of thermally activated states and concomitantly a de- 
crease in the mean free path between scattering of the 
excited states, or equivalently a decreased excitation 
lifetime.ESHSIxhus, we expect a thermally induced broad- 
ening of the excitation spectrum, as visible in Fig. [4](c). 
Below 30 K the damping is consistent with an exponen- 
tial activated behaviour T a = A a /T exp (A a /T), 
valid for gapped one-dimensional systems. 44 The fitted 
scale factors were gi = 1.05(1) and #2 = 2.22(16) for 
the low and high energy modes respectively and the fits 
are shown as solid and dashed lines in Fig. gc). A 
slight increase in excitation energy with increasing tem- 
perature was also observed [Fig. [4] (d)]. This is under- 
stood in the context of the thermally induced "blue-shift" 
that has been reported in other quasi-one-dimensional 
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systems!™l 

IV. DISCUSSION 

The excitation spectrum in BiCu2P06 is unique in 
that competing interactions along the spin chains drive 
the dynamic correlations to an incommensurate wave- 
vector. Incommensurate dynamic correlations have been 
observed in the four-leg antiferromagnetic spin tube 
Sul — Cu2ClJ^ however, the degree of incommensurabil- 
ity was much smaller than that observed in BiCu2P06- 
Furthermore, the excitations in Sul — CU2CI4 are com- 
prised of a single, triply degenerate mode, while in 
BiCu2P06 the degeneracy is split and the minimum of 
each mode occurs at a different wave vector. 

In a canonical two-leg ladder the low energy mag- 
netic excitation spectrum is composed of a well de- 
fined branch of dispersing triplets. Application of an 
external magnetic field splits the degeneracy and three 
branches of excitations may be observed.^ The degen- 
eracy may also be split in the presence of anisotropic 
interactions, such as the DM anisotropy, in the micro- 
scopic Hamiltonian!^^ We believe that the observed 
mode splitting in BiCu2P06 may be accounted for by sig- 
nificant staggered DM interactions. In BiCu2P06 local 
inversion symmetry is broken allowing a DM interaction, 
with the particular crystallographic symmetries admit- 
ting a DM anisotropy with a staggered DM vector on each 
rung as well as DM interactions on the chain bondsP^ 
Although there are many possible DM vectors, our INS 
measurements can only resolve two-bands of excitations 
in zero-field indicating that the magnetic Hamiltonian re- 
tains a local symmetry, so that we expect a single, dom- 
inant, anisotropy term. Indeed DM interactions account 
for many features of the thermodynamic measurements 
including the weak linear field dependence of the low field 
magnetization below Hj 23 l 25 l, highly anisotropic critical 
fields 25 , and absence of squar e root singularity in mag- 
netization just above H c ! 23 1 51 1 Field theoretical modeling 
has shown that staggered anisotropic interactions in a 



spin- 1/2 two- leg ladder can act to split the lowest lying 
triplet into a doublet with energy < A t and higher en- 
ergy mode with A3 > A t !^Our temperature dependent 
measurements are consistent with this scenario, although 
further detailed magnetic field dependent measurements 
are required to definitively establish the degeneracy of 
each mode. 

In summary, we have measured the spin excitation 
spectrum in BiCu2P06 at 6 K. The results confirm that 
BiCu2P06 is in a quantum disordered phase and that the 
J1-J2-J4 two-leg ladder model with interladder exchange 
J3 appropriately accounts for the dominant exchange 
pathways. The measured incommensurate wavevector 
C = 0.5747T allows for an estimate of the relative strength 
of the competing magnetic interactions J2/J1 ~ 1. How- 
ever, the lowest excitation gap is 1.90(9), is significantly 
lower than the value of 2.9 predicted by thermodynamic 
measurements. The discrepancy may be explained by 
the presence of strong anisotropic magnetic interactions 
which split the excitation spectrum into two coherent 
branches. We hope that our data will stimulate further 
theoretical work investigating the effect of Dzyloshinsky- 
Moria interactions and frustration on the spin-liquid 
ground state of a even- leg ladder. 
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